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Abstract 

The BESS-Polar spectrometer had its first successful balloon flight over Antarctica 
in December 2004. During the 8.5-day long-duration flight, almost 0.9 billion events 
were recorded and 1,520 antiprotons were detected in the energy range 0.1-4.2 GeV. 
In this paper, we report the antiproton spectrum obtained, discuss the origin of 
cosmic-ray antiprotons, and use antiproton data to probe the effect of charge-sign- 
dependent drift in the solar modulation. 
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1 Introduction 



Antiproton spectra have been measured by BESS (Balloon-borne Experiment 
with a Superconducting Spectrometer) in a series of flights from Lynn Lake, 
Canada providing reasonable statistics above 1 GeV [1,2,3,4,5,6,7]. Below 1 
GeV statistics are limited and the effects of solar modulation are greater. 
The spectrum has a distinct peak around 2 GeV, showing the characteris- 
tic feature of secondary antiprotons produced by the interaction of Galactic 
cosmic-rays with the interstellar medium. The energy spectrum of these sec- 
ondary antiprotons should decrease rapidly toward lower energies reflecting 
the kinematic constraints on antiproton production [8,9,10] and toward higher 
energies reflecting the steep power-law spectra of primary particles producing 
the antiprotons. In addition to these secondary antiprotons, there might be 
a source of primary antiprotons. Such sources have been suggested to result 
from the evaporation of primordial black holes (PBH) or from the annihilation 
of neutralino dark matter [11,12,13,14]. The spectrum of antiprotons from a 
primary source might have a peak in the energy region below 1 GeV, giving 
a flatter composite spectrum with excess flux compared to the purely sec- 
ondary spectrum [10,13]. The influence of a low energy primary peak would 
be most evident at solar minimum [12,13] and BESS (1995+1997) measure- 
ments during this period [3,4] suggested this possibility. On the other hand, 
BESS-1998 [5] and subsequent measurements [6,7], which were taken after the 
solar minimum period, are more consistent with pure secondary nature. In 
order to resolve the important questions regarding possible novel sources, we 
undertook the higher precision measurements reported here. 

In addition, one can use the charge-sign difference of antiprotons to explore 
effects of drift in the solar modulation. A detailed understanding of the effects 
of solar modulation is important to establish the existence of any primary 
sources. During the positive polarity phase of solar activity before 2000 the 
measured antiproton-to-proton (p/p) ratio [5] showed no distinctive variation. 
After the reversal of the solar magnetic field in 2000, a sudden increase of the 
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p/p ratio [6] was clearly observed. Our measurements, as well as measurements 
of the positron fraction (e + /(e + +e - )) [15], generally support recent calcula- 
tions [10,16] incorporating steady-state drift models and charge-dependent 
effects of solar modulation. Protons and antiprotons have significantly differ- 
ent interstellar spectra and the drift directions are opposite because of the 
opposite charge sign. The combination of these effects implies that the p/p 
ratio should display a more interesting evolution [16] during 2000-2010 than it 
did during the 1990's. The secondary antiprotons are also useful to probe the 
average cosmic-ray proton spectrum over a large region of the Galaxy [17]. 

The BESS-Polar experiment was proposed as an advanced BESS program 
of long-duration balloon flights over Antarctica and has been prepared since 
2001 [18,19,20] to further investigate elementary particle phenomena in the 
early Universe through a precise measurement of the low-energy antiproton 
spectrum and to search for primary antinuclei in cosmic rays. We report here 
a measurement of the cosmic-ray antiproton spectrum based on 1,520 antipro- 
ton events detected during an 8.5-day flight of the BESS-Polar spectrometer 
performed in December 2004. We discuss the origin of cosmic-ray antiprotons, 
and the solar modulation effect based on the measured p/p ratio. 



2 Spectrometer 

The BESS-Polar superconducting spectrometer, shown in Fig. 1, has been 
developed to reduce material thickness along the particle trajectory and to 
meet the severe requirements for long duration balloon flights over Antarc- 
tica [18,20]. The basic instrument concepts, such as the cylindrical configu- 
ration with an open and wide acceptance and redundant measurements for 
particle identification, are inherited from the BESS spectrometer [21]. A very 
low instrumental energy cutoff for antiprotons was achieved with a new thin- 
walled (2.46 g/cm 2 /wall including cryostat) superconducting magnet [22,23] 
using a new high-strength aluminum-stabilized superconductor. In addition, 
the outer pressure vessel was eliminated and the detectors were reconfigured. 
Low-energy particles only have to traverse 4.5 g/cm 2 of instrument material 
to be detected, about one quarter of that in the previous BESS spectrome- 
ter [21,24]. With these changes, the lowest energy for antiproton detection has 
been reduced to 0.1 GeV at the top of atmosphere (TOA). 

A uniform field of 0.8 Tesla is produced by the superconducting coil, and 
its continuous operation is maintained for over 11 days with a 400-liter liq- 
uid helium reservoir. A JET-cell type drift chamber (JET) and two inner drift 
chambers (IDCs), which were also used in the BESS-TeV spectrometer [24,25], 
are placed inside the warm bore (0.80 m in diameter and 1.4 m in length). 
The JET and IDCs are filled with C0 2 gas, and fresh gas is circulated using a 
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semi-active flow control system. Tracking of an incident particle in the x — y 
plane (perpendicular to the magnetic field) inside the JET and IDCs is per- 
formed by fitting up to 52 hit-points, each with 150 /itn resolution, resulting 
in a magnetic-rigidity (R = pc/Ze, momentum divided by electric charge) res- 
olution of 0.4 % at 1 GV, and a maximum detectable rigidity (MDR) of 240 
GV. Tracking in the z coordinate (parallel to the magnetic field direction) is 
done by fitting points inside IDCs measured by vernier pads with an accuracy 
of 1.0 mm and points in JET measured by charge division with an accu- 
racy of 37 mm. The continuous and redundant 3-dimensional tracking enables 
BESS-Polar to distinguish backgrounds such as tracks having interaction or 
scattering. A truncated mean of the integrated charges of the hit-pulses from 
the JET also provides an energy loss (dE/dx) measurement with a resolution 
of 10 %. 



The top and bottom scintillator hodoscopes, the time-of-flight (TOF) detec- 
tors, are used to determine the incident particle velocity, (3 = v/c. They also 
make two independent dE/dx measurements. The scintillator hodoscopes con- 
sist of 10 top and 12 bottom plastic scintillators with a cross section of 96.5 mm 
(width) x 10 mm (thickness) and with photo-multiplier tubes (PMTs) at- 
tached by acrylic light guides at each end. The timing resolution of each 
hodoscope is 110 ps, resulting in a f3~ l resolution of 3.3 %. In addition, a 
thin scintillator middle-TOF (MTOF) is installed on the lower surface of the 
solenoid bore to detect low energy particles which cannot penetrate the mag- 
net wall. The MTOF consists of 64 plastic scintillator bars with a cross section 
of 10 mm (width) x 5 mm (thickness) read by multi-anode PMTs. The (3 of 
the low energy antiprotons can be measured as a combination of top-middle 
TOF with a f3~ l resolution of 4.5 % below 0.6 GV (0.2 GeV for proton and an- 
tiproton). A Cerenkov counter with silica-aerogel radiator (ACC) is installed 
below the magnet. The radiator is selected to have a refractive index of 1.02 
in order to veto e~ / backgrounds up to 4.2 GeV. The top and bottom TOF 
hodoscopes and the ACC operate at ambient pressure. 



The event data acquisition sequence is initiated by a coincidence of signals in 
the top-bottom TOF or top- middle TOF hodoscopes [26]. Based on digitized 
detector information sent from the front-end electronics, the event data are 
built and recorded to onboard hard disk drives with a total capacity of 3.5 
Tera Bytes. During the Antarctic flight, the data acquisition rate was around 
1.4 kHz, and the dead time was 150 /is per event (20 % of the data-taking 
time). To supply electric power to the electronics onboard the payload, a solar 
cell array of 900 W capacity was mounted on an omni-directional octagonal 
frame around the payload. The total power consumption of the spectrometer 
was 420 W. 
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3 Data Analysis 

3.1 Balloon flight observation 

The first BESS-Polar long duration balloon flight over Antarctica was launched 
from Williams Field (77°51.8'S, 167°5.4'E) near McMurdo Station on Decem- 
ber 13, 2004 [27]. During the flight, some PMTs of the TOF hodoscopes had 
to be turned off because they were drawing excessive current, and the usable 
geometrical acceptance was consequently reduced to 73 % of the design value 
(0.157 m 2 sr at 0.2 GeV and 0.166 m 2 sr at 2.0 GeV). The flight was terminated 
on December 21 and the payload landed at the south end of the Ross Ice Shelf 
(83°6.0'S, 155°35.4'W) after a continuous observation period of 8.5 days. The 
flight trajectory was close enough to the South magnetic pole that the geo- 
magnetic cutoff rigidity was below 0.2 GV, lower than the lowest detection 
limit of the spectrometer. During the live data-taking time of 507,075 seconds 
at an average floating altitude of 38.5 km (residual atmosphere of 4.3 g/cm 2 ), 
894,482,590 cosmic-ray events were accumulated without any online event se- 
lection as 2.14 terabytes of data recorded on the hard disk drives. 



3.2 Event selection 

Since the spectrometer is cylindrically symmetric, it may be assumed that 
antiprotons behave exactly like protons in the instrument except for the sign 
of their deflection in the magnetic field and their inelastic interactions. Thus, 
all the selection criteria were defined based on the measured properties of 
protons. Initially, events were selected with (1) a single track and a downward- 
going particle fully contained in the fiducial region of the tracking volume, (2) 
only one or two hits each in the top and bottom TOF hodoscopes, (3) the hit 
position at the TOF hodoscopes consistent with the extrapolated track inside 
the JET and IDCs, (4) cosine of zenith angle of the incident particle larger 
than 0.8, and (5) either top or bottom TOF hodoscopes read by PMTs at 
both ends. As a consequence of (5), the effective geometrical acceptance had 
to be further reduced to be 40 % of the design value. 

The selection efficiency was estimated using a Monte Carlo (MC) simulation by 
applying the same selection criteria to the simulated events as to the observed 
data. The MC simulation was based on a GEANT3/GHEISHA code [28,29] 
and tuned to reproduce the results of an accelerator beam test of the previous 
BESS spectrometer [30] in which the detector configuration and materials are 
similar to the BESS-Polar spectrometer. The efficiency varied from 84.2±5.0 % 
at 0.2 GeV to 90.1±5.0 % at 2.0 GeV. The systematic error of the efficiency 
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was determined using the accelerator beam test data. Each analog-to-digital 
board for TOF has an individual dead time due to the period of switching 
reference capacitors used for baseline subtraction. Events digitized during the 
dead time have incorrect charge data and were rejected to ensure the quality 
of the data. The efficiency of surviving this cut was estimated as 77.8±0.1 % 
at 0.2 GeV and 78.8±0.1 % at 2.0 GeV. Here and in the selections described 
below, a sample of the proton data has been used to determine the efficiencies. 

The following track quality cuts were applied: (1) the reduced x 2 m the x-y 
and y-z track fitting to be less than 5, (2) the fitting error on the curvature 
(inverse rigidity) to be less than 0.015 GV -1 , (3) the track fitting path length 
to be longer than 500 mm, and (4) the residual between hit position at TOF 
obtained from the time difference in two PMTs and the extrapolated track of 
the JET to be less than 50 mm. The quality cut efficiency was estimated as 
94.7±0.2 % at 0.2 GeV and 91.7±0.1 % at 2.0 GeV. 



3.3 Particle identification 

Proton and antiproton candidates were identified with a combination of a 
charge selection and a mass selection as follows: (1) Particle charge, Z is iden- 
tified by dE/dx measurements. The dE/dx measurements with both top and 
bottom TOF hodoscopes were required to be inside a band defined as a func- 
tion of rigidity corresponding to singly-charged particles. (2) Particle mass m 
is reconstructed with rigidity R, velocity (5 and Z as m = ZeR^Jl/(3 2 - 1. It 
was required that 1/(3 be inside a band defined as a function of rigidity as 
shown in Fig. 2, so that the reconstructed mass was consistent with that of a 
proton or antiproton. The dE/dx and (3 bands for antiprotons were defined in 
the same way as for protons except for the rigidity sign. The selection efficiency 
of each dE/dx band cut was estimated as a fraction of the number of selected 
events among a proton sample selected by the other independent dE/dx se- 
lections. The net efficiency for the two dE/dx band cuts was 96.5±0.1 % at 
0.2 GeV and 96.4±0.1 % at 2.0 GeV. The selection efficiency of the 1/(3 band 
cut was estimated as 99.5±0.2 % at 0.2 GeV and 97.0±0.1 % at 2.0 GeV. 

In order to eliminate e~ and fi~ backgrounds that mimic relativistic antiproton 
candidates with (3 > 0.9, Cerenkov veto cuts were applied to select (3) the 
particle trajectory to be inside an aerogel fiducial volume, and (4) the Cerenkov 
outputs to be less than a threshold. The rejection factor was estimated as 
(8.9±0.2) xlO 2 using a relativistic proton sample with rigidity larger than 
20 GV. The efficiency including a loss of fiducial volume by 18 % was estimated 
as 75.8±0.1 % at 0.4 GeV and 62.3±0.1 % at 2.0 GeV. 

Because of the relativistic rise of dE/dx in the JET gas, ultra-relativistic 
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particles {(3 — >1) have dE/dx about 1.4 times higher than minimum ionizing 
particles. A tight cut on dE/dx measured with the JET can eliminate part of 
the e~ and yT backgrounds. (5) The dE/dx measurements with the JET were 
required to be inside a band defined as a function of rigidity. The rejection 
factor was estimated to be 1.75±0.03 at 2.0 GeV and 4.85±0.08 at 4.0 GeV, 
as a fraction of the rejected events among the e~ and \i~ sample which was 
obtained with the same selection criteria as antiprotons excluding Cerenkov 
veto cuts. The efficiency was estimated to be 98.2±0.1 % at 2.0 GeV and 
96.9±0.1 % at 4.0 GeV. 

Following the particle identification procedure, 1,512 antiproton candidates 
were identified above 0.2 GeV as shown in Fig. 2. Antiproton candidates that 
have wrongly reconstructed (3 and lie above the selection band were rejected 
by the 1/(3 band cuts. The antiproton flux was obtained by correcting for the 
cut efficiency estimated using the proton sample. The event sample obtained 
with the antiproton selection criteria excluding Cerenkov veto cuts and the 
JET dE/dx cut consists of mostly (> 99.8 %) e~ and background events. 
The number of such events divided by the total rejection factors defined above 
was used as an estimation of the number of contaminating backgrounds within 
the antiproton selection band. The estimated number of background events 
(and the ratio to antiproton candidates) was 18.7±0.3 (8.5±0.1 %) in the 1.7- 
2.1 GeV band and 8.1±0.1 (4.0±0.1 %) in the 3.4-4.2 GeV band. In the region 
-1.5 GV < R < and 0.9 < 1//3 < 1.1, the remaining background events can 
be selected with only the 1/(3 band cut and the number of such events was 
235, which was consistent with the estimated number, 243±11. 

Albedo and spillover from positive rigidity particles were negligible because of 
the high /3~ x and rigidity resolutions. To check against the "re-entrant albedo" 
background, we confirmed that the trajectories of all antiprotons could be 
traced numerically through the Earth's geomagnetic field back to the outside 
of the geomagnetic field [31,32]. 



3.4 Selection with top-middle TOF below 0.2 GeV 

Some of the antiprotons below 0.2 GeV cannot reach the bottom TOF ho- 
doscope due to energy and annihilation losses. To analyze such low energy 
events, we defined another set of selection criteria with a combination of the 
top and middle TOF hodoscopes. They are basically the same as those used 
with the top and bottom TOF, but the Cerenkov veto cuts were not applied 
since at low energies antiproton candidates are well separated from e~ and 
backgrounds using only the 1/(3 band cut. Using the top-middle TOF se- 
lection criteria, 8 antiproton candidates were identified in the energy range 
0.1-0.2 GeV. 
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Stopping antiprotons can generate secondaries in the upward direction that 
make additional tracks inside the JET. Such events are rejected by the single- 
track selection. To recover such events, a search was made for antiproton 
candidates among multi-track events. The major backgrounds for this search 
are upward moving proton secondaries generated in the lower half of the spec- 
trometer from interactions of high energy protons. Upward moving protons 
can usually be distinguished with the measured [3, but if the incident proton 
and the upward moving secondary hit the same top TOF counter, the (3 is 
not correctly measured. Such events were rejected by applying the following 
cuts: (1) Stopping antiprotons cannot generate high energy secondaries be- 
cause of the kinematical restriction. Hence for this search rigidities of all the 
tracks inside the JET were required to be less than 1 GV. (2) It was required 
that the distance between the two tracks at the top TOF to be larger than 
150 mm. The efficiency of these cuts was estimated as 72.1 ±2.1 % by applying 
the same cuts to MC generated antiproton events. A total of 31 multi-track 
antiproton candidates were found, but all of them were rejected by the cuts. 
According to the MC study, the total antiproton detection efficiency was in- 
creased by 11.1 % at 0.1 GeV by accepting multi-track events. This increase 
was treated as a systematic error, but it was smaller than the statistical error 
(+41 %/-34 %) of single-track antiproton candidates below 0.2 GeV. 



3.5 Flux determination 

After the antiproton candidates were identified, energy-dependent corrections 
were applied for backgrounds and detection efficiency. Then the absolute flux 
at the top of the instrument (TOI) was obtained by taking account of energy 
loss inside the spectrometer, live time, and geometrical acceptance. The energy 
of each particle at TOI was calculated by summing up the ionization energy 
losses inside the instrument as determined by tracing back the event trajec- 
tory. The effective geometrical acceptance was estimated using the MC code 
GEANT3 [28] and the simulation technique [33] to be 0.114±0.001 m 2 sr at 
0.2 GeV and 0.121±0.001 m 2 sr at 2.0 GeV. The error arising from uncertainty 
in the detector alignment was estimated to be 1 %. 

In order to obtain the flux at the top of the atmosphere (TOA), we applied a 
correction for survival probability of the flux reaching TOI from TOA, and a 
subtraction of the secondary component produced within the overlying atmo- 
sphere. The survival probability was estimated as 89.1±2.0 % at 0.2 GeV and 
92.4±2.0 % at 2.0 GeV based on total interaction lengths of 32.7±0.7 g/cm 2 
at 0.2 GeV and 59.9±1.2 g/cm 2 at 2.0 GeV [34,35]. 

The atmospheric secondary antiproton flux was estimated by solving simulta- 
neous transport equations [34,36]. At low energies below 1 GeV, there exists a 
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significant contribution of non-annihilating inelastic interactions or so-called 
"tertiary" antiprotons. The interaction length used was based on Stephen's 
model [35] , and the energy distribution of the tertiary production was tuned to 
reproduce the atmospheric antiproton flux measurements [37,38]. The amount 
of secondary subtraction was 11.8±1.7 % at 0.2 GeV and 28.6±4.1 % at 
2.0 GeV. The relative error of 14.3 % is composed of uncertainty in the resid- 
ual air depth (5.0 %), in the cross section of primary cosmic rays with air 
nuclei (8.9 %), and in the tertiary production (10.0 %). 



4 Results and Discussions 

We obtained the antiproton flux at the TOA in the kinetic energy range 
0.10-4.20 GeV as shown in Fig. 3 and tabulated in Table 1. The lowest en- 
ergy was determined by the detector cutoff energy, and the highest energy 
was determined by the antiproton threshold energy of the aerogel Cerenkov 
counter. The overall uncertainties including statistical and systematic errors 
are -34.1 %/+41.8 % at 0.16 GeV and ±10.9 % at 3.7 GeV with the given en- 
ergy bin width. The statistical errors are dominant over the systematic errors 
below 1.4 GeV. Owing to the long duration observation the statistical errors 
were improved from the previous measurements with BESS, carried out on 
conventional 1 or 2 day balloon flights. 

Fig. 3 also shows the results from previous BESS flights around solar mini- 
mum (95+97) [3,4] and maximum (2000) [6], which are compared with four 
theoretical calculations. The solid curves are calculations of secondary an- 
tiproton spectra with the Standard Leaky Box (SLB) model modulated with 
a steady state drift model [16], in which the modulation is characterized by 
a tilt angle of the heliospheric current sheet and the Sun's magnetic polar- 
ity (denoted as +/-). The dashed curves are calculations with the Diffusion 
plus Convection (DC) model modulated with a drift model [10,42]. Three tilt 
angles, 10°(+), 70°(— ), and 30°(— ) roughly correspond to the measurements 
with BESS(95+97), BESS(2000) and BESS-Polar(2004), respectively [43,44]. 
The dotted curves are calculations with the DC model [10] modulated with 
a standard spherically symmetric approach [39], in which the modulation is 
characterized by a single parameter (0) irrespective of the Sun's polarity. For 
each measurement, was obtained by fitting that proton spectrum measured 
by BESS, assuming the interstellar spectrum in Ref. [40]. Three values of 
0, 550 MV, 1400 MV, and 850 MV correspond to the measurements with 
BESS (95+97), BESS (2000) and BESS-Polar(2004), respectively. The dash-dot 
curves are calculations of antiproton spectra from evaporation of Primordial 
Black Holes (PBH) with an explosion rate of 0.4 x 10 _2 pc _3 yr _1 [12,41] mod- 
ulated by the same spherically symmetric approach with of 550 MV and 
850 MV. The expected signal from the PBH is modulated more than the sec- 
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ondary antiproton spectrum because of its spectral shape having the peak in 
the low energy region. While the BESS (95+97) data were suggestive of an 
excess flux at energies below 400 MeV, we do not find further evidence in the 
new BESS-Polar data presented here, even with the extended energy range. 

The p/p ratio can provide a useful probe to study solar modulation and 
its charge-sign dependence. The BESS collaboration has tracked this ratio 
through most of a solar cycle. Fig. 4 shows the p/p ratio obtained with the 
BESS-Polar measurement together with the results from previous BESS flights 
around solar minimum (95+97) [3,4] and maximum (2000) [6], which are com- 
pared with the same three calculations of secondary antiproton spectra as 
shown in Fig. 3. Fig. 5 shows time variations of p/p ratio at three differ- 
ent energies, 0.3 GeV, 1.0 GeV, and 1.9 GeV, each compared with the same 
calculations as shown in Fig. 4. 

Since the p/p ratios obtained using drift models are given as a function of tilt 
angle, they were converted into those as a function of time by taking the tilt 
angle as the mean position of the monthly variation of the maximum latitu- 
dinal extent of the current sheet [43,44]. The time variation of p/p ratios by 
the spherically symmetric approach were estimated by using a linear relation 
between 4>bess an d N C l, where 4>bess was obtained by fitting each proton 
spectrum measured with 6 BESS flights [40] and Ncl is the monthly averaged 
count rates of the Climax neutron monitor [450 




Figs. 4 and 5 show that drift models and the symmetric model reproduce 
equally well the stable p/p ratio during the positive phase. The sudden increase 
of the ratio observed by BESS measurements after the positive-to-negative 
solar field reversal is better reproduced by the drift models at energies below 
1 GeV. On the other hand, during the negative phase, where the p/p ratio 
depends on the tilt angles more than during the positive phase, the two drift 
model calculations [10,16] significantly differ from each other. We found that 
the difference mainly comes from the difference in the modulation of protons, 
which in the negative phase should come from a "horizontal" direction in 
heliosphere, i.e. along the current sheet. This would imply that the model 
of particles behavior along the current sheet is not yet established well. Our 
data prefer calculations by Bieber et al. [16] but those with the spherical 
model are also preferred. In addition, Fig. 3 shows that the spherical model 
reproduces antiproton spectra of all the three BESS measurements better than 
drift models. This implies that the charge-sign dependence in the modulation 
is not so significant during the negative phase. The spherical model is still the 
best one to reproduce measurements. 



U.S. National Science Foundation Grant ATM-0339527 




10 



5 Conclusion 



The first BESS-Polar experiment was carried out in Antarctica in December 
2004, a transient period before the solar minimum in 2007. Using a new spec- 
trometer with reduced material thickness and a long duration balloon flight 
near the Earth's south magnetic pole, the lowest energy limit of the antiproton 
flux was extended down to 0.1 GeV, and the statistics were improved com- 
pared with the previous BESS experiments. The series of BESS measurements 
have enabled a crucial test of models of the solar modulation. Drift models can 
reproduce the drastic behavior of the p/p ratio around the solar field reversal, 
but for the negative phase the spherical model is still best able to reproduce 
most of the measurements. BESS data should motivate further development 
of the drift models with more realistic parameters and time dependence. For 
a more extensive future search for cosmic-ray antiprotons of primordial ori- 
gin the BESS-Polar result provides an important baseline measurement of the 
secondary antiproton spectrum, which will be compared with a spectrum mea- 
sured in December 2007 through January 2008, by the second Antarctic flight, 
during the solar minimum period. 
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Table 1 

Antiproton flux and p/p ratio at the top of atmosphere with statistical (first) and 
systematic (second) errors. Np and Nbg are the number of observed antiprotons 
and estimated background events, respectively. 



Kinetic energy (GeV) 
range mean 



iVp N BG 



p flux 

( m -2 sr -l s -l Ge y-l) 



p/p ratio 



0.10-0.18 
0.18-0.28 
0.28-0.40 
0.40-0.56 
0.56-0.70 
0.70-0.88 
0.88-1.10 
1.10-1.37 
1.37-1.72 
1.72-2.15 
2.15-2.68 
2.68-3.36 
3.36-4.20 



0.16 
0.22 
0.34 
0.48 
0.65 
0.79 
1.00 
1.23 
1.54 
1.92 
2.40 
3.01 
3.68 



15 

24 
23 
37 
55 
84 



276 
204 



0.0 
0.0 
0.0 
0.1 
0.2 
0.6 
3.3 



143 10.8 

198 18.2 

220 18.7 

233 15.1 



11.9 
8.1 



4.37 



+1.81+0.25 
-1.47-0.25 



x 10- 



4-53ll:Slg:i? x 10-3 
6-67±i:«±8:i? x 10-3 

X 10-3 



o Q/1 +0.91+0.30 
°- y ^-0.80-0.30 



8.27 



+1.53+0.47 
-1.39-0.47 



Q 8O +1.41+0.59 
y - oz -l. 29-0.59 

1 17 +0.l4+0.07 
L - L '-0.13-0.07 

+0.15+0.12 
-0.14-0.12 



1.64 



-, o 4 +0.15+0.13 
i -°^-0.14-0.13 



1.62 



+0.11+0.13 
-0.11-0.13 



1 /m+0.09+0.12 
1 - flu -0.09-0.12 



1.47 



+0.09+0.11 
-0.09-0.11 



1 i n +0.08+0.09 
i - lu -0.08-0.09 



X 10-3 
X 10-3 

x 10- 2 
x 10- 2 

io- 2 
io- 2 
io- 2 

x IO" 2 
x IO" 2 



X 



X 



8.28+| 
9.49l^ 
5.181}; 

i-iolg: 

1-3818; 

1.831°; 
2.94l° 
4.0018; 
4.4218 
5.0118; 
7.1718; 
7.4618; 



03+0.64 
39-0.64 

85+0.40 
10-0.40 

06+0.38 
20-0.38 

18+0.06 
20-0.06 

18+0.08 
20-0.08 

21+0.11 
22-0.11 

26+0.20 
27-0.20 

31+0.28 
33-0.28 

30+0.33 
30-0.33 

33+0.41 
33-0.41 

43+0.54 
43-0.54 

52+0.60 
52-0.60 



X IO" 6 
X IO" 6 
X IO" 6 

x IO" 5 
x IO" 5 
x IO" 5 
x IO" 5 
x IO" 5 
x IO" 5 
x IO" 5 
x IO" 5 
x IO" 5 
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Incoming particle 




Fig. 1. Cross-sectional view of the BESS-Polar spectrometer. 
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Fig. 2. The versus rigidity plot, and antiproton selection band. The same band 
but opposite rigidity sign is applied to select protons. For the negative rigidity, 
all the events after Cerenkov veto cuts and JET dE/dx cut are shown. For the 
positive rigidity, 0.5 % of the events after Cerenkov veto cuts are shown. Antiproton 
candidates are shown with dots of larger size. 
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10 1 1 

Kinetic energy (GeV) 

Fig. 3. Antiproton flux at the top of the atmosphere obtained with the first BESS-Po- 
lar flight together with results from previous BESS flights around solar minimum 
(95+97) [3,4] and maximum (2000) [6]. The solid curves are calculations of sec- 
ondary antiproton spectra with the Standard Leaky Box (SLB) model modulated 
with a steady state drift model [16] by solar tilt angles and magnetic polarities of 
(from top to bottom, the first two are very close) 10°(+), 10°(-), and 70°(-). The 
dashed curves are calculations with the Diffusion plus Convection (DC) model [10] 
modulated by (from top to bottom, the first two are very close) 10°(+), 30°( — ), 
and 70°(— ). The dotted curves are calculations with the DC model [10] modulated 
with a spherically symmetric model [39] by (from top to bottom) 550 MV, 850 MV, 
and 1400 MV. The dash-dot curves are calculations of antiproton spectra from 
evaporation of primordial black holes with an explosion rate of 0.4 x 10~ 2 pc~ 3 yr _1 
modulated by 550 MV(top) and 850 MV(bottom) [12,41]. 
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Fig. 4. p/p ratio obtained with the first BESS-Polar flight together with results 
from previous BESS flights around solar minimum (95+97) [3,4] and maximum 
(2000) [6]. The solid curves are calculations with the Standard Leaky Box (SLB) 
model modulated with a steady state drift model [16] by solar tilt angles and mag- 
netic polarities of (from bottom to top, the first two are very close) 10°(+), 10° (— ), 
and 70°(— ). The dashed curves are calculations with the Diffusion plus Convection 
(DC) model [10] modulated by (from bottom to top, the last two are very close) 
10°(+), 30°(-), and 70°(-). 
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Fig. 5. The top panel shows time variations of the tilt angle [43,44] shown by a 
solid curve and modulation parameter for the spherically symmetric model [39] 
shown by a dotted curve. The modulation parameter values continuous over time 
were estimated by using its linear relation with Climax neutron monitor data [45]. 
The linear relation was established using the modulation parameter for each BESS 
flight obtained by fitting the BESS proton spectrum. The other three panels show 
time variations of the p/p ratio at 0.3 GeV (2nd), 1.0 GeV (3rd), and 1.9 GeV 
(bottom). The data of p/p ratio are compared with time variations predicted by 
two drift models shown by solid curves from Bieber et al. [16] and dashed curves from 
Moskalenko et al. [10], and with the spherically symmetric moduation by Fisk [39] 
shown by dotted curves. 
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